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We present a theory for the memory effect in electron glasses. In fast gate voltage sweeps it is 
manifested as a dip in the conductivity around the equilibration gate voltage. We show that this 
feature, also known as anomalous field effect, arises from the long-time persistence of correlations in 
the electronic configuration. We argue that the gate voltage at which the memory dip saturates is 
related to an instability caused by the injection of a critical number of excess carriers. This saturation 
threshold naturally increases with temperature. On the other hand, we argue that the gate voltage 
beyond which memory is erased, is temperature independent. Using standard percolation arguments, 
we calculate the anomalous field effect as a function of gate voltage, temperature, carrier density 
and disorder. Our results are consistent with experiments, and in particular, they reproduce the 
observed scaling of the width of the memory dip with various parameters. 

PACS numbers: 73.61. Jc 



I. INTRODUCTION 

In the insulating low temperature phase of dirty semi- 
conductors or granular metallic films the unscreened 
Coulomb interactions between localized electrons lead to 
glassy behavior such as slow relaxation, history depen- 
dence of observables, non-ergodicity, and memory effects. 
Even though such "Coulomb glasses" were theoretically 



predicted more than twenty years age 



.1.2.3.4 



', it was a ma- 



jor task to provide convincing experimental evidence for 
their existence. To our knowledge, very slow electronic 
relaxation was first reported in the context of capaci- 
tance measurements in doped GaAs by Monroe et al£. 
At temperatures well below IK, they observed relaxation 
times that reached the scale of seconds. Even more strik- 
ing non-equilibrium behavior was found by Ovadyahu's 
group in the conductivity of strongly disordered indium- 
oxide films&, where the logarithmic relaxation can ex- 
tend up to several hours (a typical experimental set-up 
is shown in Fig. Over the last decade, careful stud- 
ies of these systems have demonstrated that the elec- 
tronic out-of-equilibrium behavior is indeed due to the 
strong frustration induced by the Coulomb interactions 
between localized electrons, and does not primarily re- 
flect the glassy dynamics of extrinsic degrees of freedom 7 . 
All the key features usually associated with a glassy sys- 
tem have been observed in these systems: Aging 8 ^, the 
dependence of sample properties on its historjiiSiii and 
memory effects^. The latter appear as a dip in the film 
conductivity as the gate voltage is swept through the 
point at which the system was equilibrated for a long 
time. The memory of these equilibration conditions usu- 
ally persists for several hours after the gate voltage has 
been changed to a new value. 

A very similar anomalous field effect, accompanied by 
slow relaxation, was observed in various granular metals 
such as Aui£, Ali^i^, as well as Bi and PblSJI. Fur- 
thermore, the aging behavior and the temperature de- 
pendence of the memory dip reported in granular A! 18 
are very similar to those found in indium-oxide films 19 . 
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FIG. 1: Typical experimental setup: The sample is a film 
(semiconductor or granular metal) of thickness d, which is 
coupled capacitively to a gate electrode. Variation of the gate 
voltage V g slightly changes the number of carriers in the film. 
The conductivity of the sample is probed through contact 
electrodes that are directly attached to the sample. 



This suggests that these glassy effects are rather univer- 
sal, even though the details of the hopping mechanism 
and the temperature dependence of the resistivity are 
clearly different in the two systems. 

The anomalous resistance fluctuations observed in ul- 
trathin granular aluminium nlmsi^ and silicon MOS- 
FETs close to the metal insulator transition22i2i were 
also interpreted as indications for glassy behavior. Un- 
fortunately, in these systems, it is difficult to disentangle 
effects due to intrinsically glassy behavior of interact- 
ing electrons from the strong response of the percolating 
network of hopping electrons to extrinsic slow degrees of 
freedom. 

It has been conjectured that the glassy memory dip 
reflects a Coulomb gapSS in the density of states which 
arises from the unscreened Coulomb interactions. This 
conjecture is rather natural in the light of recent theoret- 
ical work which suggests that the emergence and the uni- 
versality of the Coulomb gap are related to a glass tran- 
sition at a finite T g in these systemsSSiSSiSl. This opens 
the appealing perspective of obtaining more detailed in- 
formation on Coulomb correlations from a quantitative 
analysis of the memory experiments. 



So far, the precise connection between the memory ef- 
fect and Coulomb correlations in the electron configu- 
ration has remained unclear. The aim of this paper is 
to provide a quantitative analysis of this effect, assum- 
ing a number of glassy properties of the electron system. 
We start with a review of the glassy features observed in 
experiments (Section [nj , concentrating on the memory 
effect. In Section IlIII we briefly summarize the theoreti- 
cal background on the Coulomb gap, hopping transport, 
percolation theory and glassy behavior which is needed 
for the quantitative theory of the memory effect in Sec- 
tion The approximations underlying our theory and 
some open issues are discussed in Section [V] We con- 
clude with a brief summary of the main results. Several 
detailed discussions have been deferred to appendices in 
order not to interrupt the main line of the reasoning. 



II. GLASSY BEHAVIOR IN EXPERIMENTS 

A. Slow relaxation and aging 

In this section, we discuss some of the experiments ex- 
hibiting glassy behavior. Monroe et al£ were the first to 
observe slow electronic relaxation in the capacitance of p- 
type doped, partially compensated GaAs. More recently, 
very slow logarithmic relaxation was observed in the con- 
ductivity of various granular metalsi£*i£ii£, as well as in 
indium-oxide films 7 . Furthermore, after equilibration un- 
der fixed experimental conditions and subsequent moder- 
ate excitation during a time t w (e.g., by gate voltage^ or 
electric field^), the typical relaxation time of such films 
is found to scale with t w . This phenomenon, known as 
simple aging, is observed in many glassy systems such 
as polymers and supercooled liquids^, as well as spin 
glasses^. 
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FIG. 2: Left: In an electron glass, the conductivity exhibits a 
symmetric dip (on top of the linear normal field effect) around 
the equilibration gate voltage Vg. Far from V g ° the conduc- 
tivity is higher since a fast change of gate voltage takes the 
system into a higher-lying metastable state. Below a charac- 
teristic temperature T g , there is a clear difference between the 
conductivity in an "equilibrated" state (EQ) and metastable 
states (NEQ), see the inset. The relaxation from a metastable 
state to an "equilibrated" state requires collective electron 
hops (lower left pannel) or crossing of high thermal barriers. 
Both are so slow that they cannot entirely take place on ex- 
perimental time scales. However, partial relaxation results in 
the decrease of the dip amplitude with the sweep rate. 
Right: An analogous experiment in spin glasses measures the 
magnetization as a function of magnetic field H after a slow 
quench in the equilibration field H°. These traces exhibit a 
peak around -ff , reflecting that the field cooled magnetiza- 
tion is larger than the magnetization that one obtains after a 
quench in H ^ H° and a subsequent switch to H°. The inset 
illustrates this point with the well-known difference between 
field cooled (FC) and zero-field cooled (ZFC) magnetization 
as a function of temperature. The relaxation from the ZFC 
state to the FC state requires the collective rearrangement of 
many spins (lower right pannel). 



B. The memory effect 

One of the most striking manifestations of the electron 
glass is the anomalous field effect: After the equilibra- 
tion of the sample at some gate voltage V®, subsequent 
traces of conductivity as a function of gate voltage keep 
a long-lasting memory of these equilibrium conditions in 
the form of a symmetric dip around V®. This dip is 
superposed on the linear normal field effect due to the 
increase of carriers, which is usually subtracted and will 
not be considered further here. In Fig. we illustrate 
the memory effect - and for the reader more familiar with 
spin glasses, it is juxtaposed with an analogous experi- 
ment that could be done in spin systems. 

The fact that the conductivity increases no matter 
whether carriers are added or depleted, can be under- 
stood on a qualitative level by the observation that any 
perturbation taking the system out of equilibrium must 
lead to an increase of the conductivity^. On a more quan- 
titative level we will have to explain the following exper- 



imental observations extracted from extensive studies on 
indium-oxide filmsiS: 

(i) The width T of the memory dip (measured as a 
function of density of induced carriers) is remarkably uni- 
versal. In particular, it is independent of the sweep rate 
or the application of a magnetic field. Even more surpris- 
ingly, it remains unchanged under thermal annealing, a 
process which reduces the disorder significantly and thus 
increases the conductivity by several orders of magnitude. 

(ii) r increases with carrier density. 

(Hi) r increases roughly linearly with temperature, and 
keeps a certain memory of temperature: After a sudden 
quench from the equilibration temperature T to a lower 
temperature T' < T, the width T relaxes only slowly to 
the value corresponding to T", keeping a memory of the 
larger width characteristic of T— . 

Note that (ii) is a strong indication for the relevance 
of electron-electron interactions!. Similar results, in par- 
ticular the decrease of T with temperature and memory 
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of higher temperatures, were recently found in granular 
aluminium films, too 1 ^ 1 ^. 

We will show in Section Hvl that these key features can 
be understood semi-quantitatively with relatively simple 
arguments on the glassy free energy landscape and the 
stability of its local valleys. 



which maps the problem onto the effective model QJ, 
as discussed in App.[H] We expect that the lattice model 
remains a reasonable approximation up to z « 1. Typi- 
cal values for z in indium-oxides can be estimated to be 
of the order of z » 0.2 — 0.5. 



III. THEORETICAL INGREDIENTS 
A. Model 

In Anderson insulators, the unscreened Coulomb in- 
teractions between the localized carriers are the crucial 
ingredient which lead to strong electron-electron corre- 
lations, the formation of the Coulomb gap and glassy 
behavior. An approximate description of such systems is 
given by the classical lattice Hamiltonian 23,24 



e riiUj 



(1) 



where tii = 0, 1 are the occupation numbers of randomly 
positioned lattice sites i, and we fix the average occupa- 
tion to 1/2. k denotes the host dielectric constant. The 
disorder energies are considered as independently and 
identically distributed random variables with a charac- 
teristic width W. This corresponds to a "bare" density 



of states 



l/r~W where 



-1/3 

- n t 



is the mean 



interparticle spacing in a system with carrier density n c . 

Two standard models for semiconductors with local- 
ized electrons are usually considered in the literatur o 31 ' 32 : 
the "classical impurity band" (CIB) and amorphous semi- 
conductors (AS). They predominantly differ in the bare 
density of states vq, see App.0 In the CIB the disorder 
W is due to randomly distributed charged impurities, so 
that the disorder is of the order of the nearest neighbor 
interactions between carriers, and 7 = W/(e 2 / nr c ) ~ 1. 
In the AS, the disorder is due to strong local inhomo- 
geneities, and can be much larger than the nearest neigh- 
bor interaction, 7> 1. In indium-oxides, one can tune 
from the CIB-regime (low n c ) to the AS-regime (high 
n c ) by controlling the carrier density n c , the crossover 
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occurring around nx 

Most interesting glassy effects have been observed in 
systems with a high density of carriers: in various granu- 
lar metals and, most importantly, in indium-oxide. The 
latter is a highly disordered semiconductor, that can be 
prepared to have exceptionally high carrier densities in 



a large range n c 



10 



10 



10 r cm , while still being 



insulating. It is not obvious that these high density sys- 
tems can still be described by the lattice model Q , since 
typically the number of carriers per localization volume 
is larger than one. However, provided the parameter 
z = v £ i D ~ 1 e 1 1 k is small (£ being the localization length), 
we can consistently restrict ourselves to a subsystem of 
well-separated electrons sitting on sites with small e,, 



B. Coulomb gap and hopping transport 

Since the pioneering works by Pollak 22 , Efros and 
Shklovski: 2 - in the early seventies, it has been known that 
the unscreened Coulomb interaction in Anderson insula- 
tors lead to important correlations in the configuration 
of electrons, and in particular to the Coulomb gap in the 
density of states. An upper bound for the single parti- 
cle density of states in D dimensions is obtained from a 
self-consistent stability argument 23-24 
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where 



E c = ({e 2 /K) is /a D )°- 



E<E C (2) 
E > E c ' 



(3) 



is the typical scale below which Coulomb correlations 
dominate over the disorder, and is a numerical con- 
stant. In J2Jl, Ei = dH/drii is the energy cost to change 
the occupation on the site i. A very similar distribution 
of energy costs (but with a substantially larger dflS) 
holds for the quasiparticle excitations relevant for vari- 
able range hopping, sometimes referred to as electronic 
polarons^ii^ 3 .. At higher temperature, the Coulomb gap 
fills in gradually, and is essentially smeared out for T > 
Ec , even though a small depression due to Coulomb cor- 
relations should persist. 

The presence of a Coulomb gap at low temperatures 
leads to a crossover of the variable range hopping con- 
ductivity form Motts's law 




with T 



(D) 
M 



R(T) = R exp 



1 / vq^ d to the Efros-Shklovskii law 



(4) 



R(T) = R exp 




(5) 



with TgJ ~ e 2 / The prefactor Rq is a slowly varying 
function of temperature. 

One can obtain a quantitative description of variable 
range hopping from the standard application of perco- 
lation theory to a network of Miller- Abrahams resistors 
formed by pairs of sites^^^SiSi, as reviewed in App. [D] 
In particular, considering an equilibrium quasiparticle 
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density of states of the form J5J one finds a crossover 
function 



log[R(T)/R ] = z- l ^ D -^Tl{T/T x ), 



where 



Tx = {tE^/T'm 



+ 1 /T?, s ) 1 /(£ , -l) 



A/{D-i) E( 



(6) 



(7) 



is the temperature where the conductivity crosses over 
from Mott's regime to the Efros-Shklovskii regime (see 
j^ e f g| 38 i 39 f or s i m ii ar approaches). In Section llVl we will 
be concerned with the modification of the conductivity 
as the density of states is driven out of equilibrium by 
the application of a gate voltage. 

For a long time, experimental evidence for the 
Coulomb gap in doped semiconductors was only indirect 
in the form of the Efros-Shklovskii hopping law, and it 
is difficult to extract detailed information on Coulomb 
correlations from the temperature dependence of the re- 
sistivity R(T) alone. However, differential measurements 
represent direct fingerprints of correlations since they are 
only sensitive to changes in the electron configuration as 
an external parameter is varied. 

In the last ten years, several tunneling experiments on 
weakly insulating samples provided direct evidence for 
a pseudogap in the density of states around the Fermi 
level^2i£ii22i2iM!4£i2§ii£. However, such experiments are 
restricted to the regime relatively close to the metal- 
insulator transition. From this point of view, the analy- 
sis of the anomalous field effect represents a convenient 
method to probe Coulomb correlations also deeper in the 
insulating regime. 



C. Glassiness 

A key ingredient to the understanding of the mem- 
ory effect is the glassiness of the electrons at low tem- 
peratures. Such a behavior can be expected from the 
theoretical consideration that Coulomb systems are very 
similar to frustrated antiferromagnets, for which recent 
experiments have demonstrated the existence of a ther- 
modynamic phase transition and spin-glass-like out of 
equilibrium behavior in almost pure samples 48 . 

From numerical simulations, it is well established that 
at sufficiently low temperatures there is a multitude of 
metastable states, which are not ergodically connected 
within timescales accessible in a simulation, since they 
arc separated by large activation or tunneling barriers 
49 i 50 i 51 i 52 i 53 i 54 . Extrapolating to experimentally relevant 
timescales, one expects that a dynamical glass transition 
takes place in real systems as well. 

More theoretical insight can be gained from mean-field 
theory^SiSi. Let us first discuss the case of strong disor- 
der (7 > 1), which corresponds to most experiments in 
indium-oxide (n c > nx)- In 3D systems, a locator ap- 
proximation to the high temperature expansion predicts 



a glass transition at finite temperature^ 

1 



T (3D) = 

9 6(2/tt) 1 /4"3 



V 2 c 
a, E c , 



(8) 



valid for large disorder, 7 > 1. Applying the same ap- 
proach to 2D systems in the strong disorder limit, one 
obtains the prediction 



T (2D) 

g 



log 



a 2 E^ D \ 



(9) 



It is possible that in systems in low dimensions the 
sharp mean field thermodynamic transition is rounded 
due to activation over finite but high barriers. In this 
event, T g in Eqs. (|8I9(I is expected to mark the crossover 
to strongly activated dynamics. 

Apart from predicting T g , mean field theory tells us 
that in 2D the ratio T g /Ec © can be numerically large, 
in particular, if we remember that the relevant value of 
cti is the one associated with the quasiparticle density 
of states. One may therefore expect glassy behavior in 
strongly disordered films in a substantial range of tem- 
peratures T/Ec > 1 where the Coulomb gap is hardly 
developed yet. It is indeed not unusual that the glass 
transition occurs at a temperature where the density of 
states does not yet exhibit any of its prominent low tem- 
perature features. This is for example the case in the 
Sherrington-Kirkpatrick model of long range spin glasses 
where a linear pseudogap starts to open only well below 
T g . The same is predicted by the mean field solution for 
strongly disordered Coulomb glasses. 

In the case of moderate disorder, 7 « 1 (the CIB 
model) mean field theory predicts a rather low glass tran- 
sition temperature 27 in 3D, 



rp(3D) 

9 



0.03 



2 1/3 



(10) 



consistent with the small values found in simulations on 
irregular lattices without on-site disordes££*£&. We expect 
a similar situation in weakly disordered 2D systems. 

In indium-oxide films, a rapid drop of the relaxation 
time is observed as the carrier density is decreased be- 
low n cr 10 19 cto~ 3 , while keeping the temperature con- 
stants 7 .. It is possible that this is a manifestation of the 
glass transition. Indeed, such samples are in the classical 
impurity band regime (n cr < nx), and Eq. (|10f> yields 
a value of T g close to the measurement temperature of 
T m = AK. Films with even lower density, n c < n cr , are 
still in their ergodic high temperature phase at T = T m , 
and the relaxation times are unmeasurably fast. 



IV. THEORY OF THE MEMORY DIP 

A. The density of states as a function of T and V g 

Our theoretical approach is based on the assumption 
that the metastable states visited after an excitation 



FIG. 3: A sketch of the density of states in a 2D electron glass. 
The dashed curve corresponds to zero temperature. The dot- 
ted curve shows the result of thermal smearing, cf. Eq. ll7> . 
The solid curve represents the density of states immediately 
after applying a gate voltage, cf. Eq. l20ll . A/i is the shift of 
the chemical potential due to the charging of the sample. 



by gate voltage reflect the way in which this state was 
reached. In particular, we argue that at a new gate 
voltage V g the (quasiparticle) density of states, and thus 
the hopping conductivity, will be distinct from the equi- 
librated state at the same V g , even if all spontaneous 
single-particle relaxations had time to take place. The 
full relaxation to equilibrium will involve multi-particle 
relaxations and/or processes with high activation ener- 
gies. Analytical and numerical arguments in support of 
this scenario have been discussed in Refi££. 

Here we are focusing on truly 2D systems, that is, 
films of thickness d < d cr = [v^e 2 / k) -1 / 2 (see Add. IU1 
for a discussion of the crossover to 3D systems). As 
discussed in Section IIII Bl the density of states in such 
films exhibits a linear Coulomb gap at low temperatures. 
At energies larger than the Coulomb correlation scale 

Eq = (e 2 / ' n) 2 v^ D ^ /a 2 , the density of states approaches 

the constant bare density of states v^ ^ = vod. Wc 
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describe this crossover (for T 
function 



0) by the interpolating 



Po(E) 



ta,nh\E/E c 



(11) 



whose precise form is, however, not essential for the fol- 
lowing analysis. 

At finite temperature some electrons are excited out 
of their local equilibrium position, which induces fluc- 
tuations in the (single electron) site energies, Ei = 



(0) 



where 



E r = E. 



(12) 



(13) 



T/E r 



FIG. 4: (Color online) The density of thermally excited carri- 
ers, ut, as a function of temperature, or, equivalently, the gate 
voltage difference CVins/e at which the anomalous held effect 
saturates. At high temperature, Vi ns is linear in T [oc (T— To)], 
and scales as T 2 at low temperatures, nr and CVi w /e are 
plotted in units of zEcv^d. 



is the energy cost to change the occupation of site i in 
the locally stable configuration characterized by {n\ }, 
and 



^e 2 8n J /nr l 



(14) 



are potential fluctuations due to thermally activated 
changes in the occupation Srij. In a first approxima- 
tion we assume the <50j's to be independent Gaussian dis- 
tributed variables with variance (5(f) 2 } = aTriT(e 2 / k) 2 , 



Pt{84>) 



1 



e 2, TjirotTriT 



exp 



OiTTlT 



(15) 



where ay = 0(1) is a numerical factor and is the 
density of thermally excited electrons 



nr= I dE p (E) (1 - f(E)) + / dE p (E)f{E), 

(16) 

/ being the Fermi distribution. Note that at high tem- 
peratures ut is linear in T, while at low temperatures it 
approaches zero as nr ~ T 2 (see Fig. 0J). 

As a consequence of these fluctuations, the density of 
states is smeared. In the approximation of independent 
shifts 8 <pi, it is described by the convolution 



p(E,T) = / d(6<j>)Pr(6<j>)po{E 



(17) 



Upon application of a gate voltage V g , new carriers are 
introduced into the system. It is reasonable to assume 
them to rapidly occupy the empty sites with the lowest 
energies Ei available. Sometimes, the occupation of a 
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site may cause a small number of neighboring particles 
to hop slightly away, to reduce the energy of the system. 
Considering the introduction of the new particle and such 
local rearrangements as one composite process, one may 
say that the new carriers actually occupy "quasiparticle" 
states (very similar to the electronic polarons relevant for 
conductivity^). However, even though the thus reached 
state may be stable to single particle relaxations, it will 
in general be an excited state under the new gate voltage. 
Only if the system is given a long time to equilibrate, it 
will relax to the new ground state, which involves mul- 
tiparticle transitions or the crossing of high activation 
energies. 

In order to describe the quasiparticle density of states 
analytically, we assume that the latter type of relaxation 
processes has not had time to occur. Furthermore, we 
assume that apart from the local response of nearby par- 
ticles, the introduction of new carriers does not trigger 
major rearrangements of the electron configuration. This 
will be justified further below. A more thorough investi- 
gation of these assumptions can be found in££. 

Under these assumptions, the effects of the gate volt- 
age are twofold: (?) the new carriers successively fill the 
empty (quasiparticle) levels close to the Fermi energy, 
shifting the chemical potential to /i 4- A/x, while the min- 
imum in the density of states remains at the old value of 
[i (cf. Fig. I2J); (ii) the extra particles further smear the 
density of states, similarly to the thermal effect described 
above. They induce further energy shifts which we 
take to be randomly distributed according to 



Pv g {54>) = 



1 



e 2 ^-KotvCVg/e 



exp 



k 2 (8cj> - Afi) 2 
"e 1 a v CVg/e 



(18) 

where ay = O(l) and C is the capacitance per unit area. 
We have also accounted for the global shift in chemical 
potential A/j,, which is related to the gate voltage by 



CV g = e 



A M (V S ) 



dEp(E). 



(19) 



Notice that in the presence of a Coulomb gap, the depen- 
dence of A/i on V g is non-linear. The density of states 
after a sudden gate voltage change is finally obtained as 



p{E,T,V g \ 



d{8<t>)PvA5cj ) )p{E-5(l>,T,V g =Q). (20) 



The density of states at different stages of smearing is 
shown in Fig EI Below, we will use the density of states 
(|20|l to calculate quantitatively the out of equilibrium 
conductivity 

Note that in assumption (i) it is implied that new car- 
riers will occupy sites across the whole film. This is only 
justified if the film thickness is of the order of the screen- 
ing length of the sample. In the absence of a Coulomb gap 
the latter can be estimated as l sc ~ (v^e 2 / k) -1 / 2 which is 
of the same order as the thickness d cr which governs the 
crossover to a 3D system. (However, in the presence of 



an Efros-Shklovskii Coulomb gap, the screening length is 
probably substantially larger on intermediate timescales, 
as suggested by the capacitance experiments of Ref£.) 
We believe that in films thicker than l sc , the anomalous 
field effect is mostly due to the filling of states within a 
screening length from the surface. The quantitative the- 
ory below does not strictly apply to this case. However, 
the instability argument given in the following subsection 
should still hold, provided the film thickness d is replaced 
by d cr ~ he- 



B. Instability criterion and breakdown of memory 

The description l|20|) of an adiabatic response to the 
change of gate voltage V g , without any relaxation of the 
electron configuration, is applicable only for small enough 
values of V g . As the gate voltage is increased, more 
and more new particles are introduced into the sample 
and reshuffle the site energies, until at a certain scale 
(Vg = Vi ns ) the local minimum in which the system re- 
sides becomes unstable. For higher gate voltages the sys- 
tem will relax to a new local minimum whose density 
of states is no longer described adequately by the adia- 
batic smearing and shifting of Eq. I|20l) alone. As long as 
the Coulomb gap is not strongly developed, it is reason- 
able to assume that the new local minimum represents a 
rather generic metastable state relatively high up in the 
energy spectrum of all possible states. We then expect 
that the conductivity will not significantly change upon 
further increase of the gate voltage, since the system re- 
mains in the high energy spectrum of states. However, 
such a new metastable state will still have a large config- 
urational similarity (or "overlap" in spin glass language) 
with the original ground state: Most of the sites that 
were occupied in the original local state remain occupied 
in the new local minimum. The memory of the original 
configuration is thus preserved. In particular, when the 
gate voltage is swept back to its original value, the low 
equilibrium conductivity will be recovered. As the gate 
voltage is increased beyond the instability scale v\ ns , the 
overlap of the new local minimum with the original state 
continuously decreases and the memory of the original 
state is gradually lost. This will be manifested by the 
disappearance of the memory dip once the gate voltage 
is swept beyond a scale V mem . 

Let us now analyze in more detail what determines 
the instability scale V5 ns - It is reasonable to expect that 
as long as the density of carriers introduced by the gate 
is smaller than the density of thermally excited electrons 
tit, the gate voltage effect is perturbative, which justifies 
our adiabatic treatment of the density of states. This 
reasoning implies Cv\ ns ^ eriT- 

On the other hand, at gate voltages CV g > eriT the 
shift in chemical potential is of the order of the tempera- 
ture Ap, ~ T, and accordingly, the new carriers are intro- 
duced on sites that were essentially empty in the original 
state. The local environment of those sites will gener- 
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FIG. 5: (Color online) Memory dip as a function of gate volt- 
age for z — 0.4 and various temperatures of the order of 
Ec and higher. We plot the relative change in conductiv- 
ity with respect to its asymptotic value at large gate volt- 
ages, AG/G = [G(V g ) - G(V inB )]/G(V ins ). The cusp width is 
proportional to temperature, and its amplitude decreases as 
temperature is increased. CV g is plotted in units of ezEcvod 



FIG. 6: (Color online) AG/G = [G(V g ) - G(V inB )}/G(Vi nB ) 
plotted against gate voltage for low temperatures, T < Ec, 
and z — 0.4. In this regime, the cusp width increases quadrat- 
ically with temperature. The adiabatic percolation treatment 
predicts the amplitude of the cusp to increase with tempera- 
ture. This is probably an artifact, see the discussion in Sec- 
tion CV g is plotted in units of ezEcvod. 



ally not be favorable to the addition of a new particle. 
Rather, the newly introduced electron will trigger fast 
relaxation processes and destabilize the original state. In 
other words, the configuration generated by occupying 
more and more levels will soon become a generic high 
energy state for CV g > enx- In summary, we expect an 
instability and thus a saturation of the out-of-equilibrium 
conductivity at 



Vn; 



, ( 7T 2 T 2 /6E, 



C, 



c 



c c 



2\n2T-n/2is d, T>E C 



(21) 



where = J dE(i>od — p(E)) ~ vodEc is the total deficit 
of density of quasiparticle states due to the presence of 
a Coulomb gap. In Fig. 0] we plot V; ns as a function of 
temperature. 

As explained above, we expect the memory dip to sat- 
urate around the instability scale, its total width being 
roughly T = 2V\ ns . A linear high temperature behav- 
ior r ~ T — Tq was indeed observed experimentally 58 . 
Equation l|21|) further predicts the interesting relation 
T = J dE(l - p(E)/v Q d)/4\n2, which allows to deter- 
mine experimentally the width of the Coulomb gap. 

The memory of the equilibrium state is essentially 
erased once the gate voltage has exceeded a certain scale 
V mcm > Vins- We expect this crossover to occur when the 
random energy shifts 8<pi are comparable to the Coulomb 
correlation scale, (Sip 2 ) 1 / 2 « Ec- More explicitly, we ob- 
tain the estimate 



V m 



Qv—^-Ec, 



(22) 



which is temperature independent contrary to Vi ns . The 
numerical factor £y = 0(1) is presumably numerically 



large (at least of the order of T g /Ec which may be ap- 
preciable in strongly disordered 2D systems, cf. (JSJ). 

In Refiii the authors reported that the ratio V mein /Vi nB 
is not a universal number (at fixed temperature), but in- 
creases with carrier density. The above arguments indeed 
suggest that at high temperatures, T > Ec, 



V m 



oc 



E c 
T-T 



(23) 



At fixed temperature, this ratio increases with carrier 
concentration, as Ec does. 



C. Non-equilibrium conductivity 

For gate voltages V < Vi ns we can calculate the non- 
equilibrium conductivity from the modified density of 
states, Eq. I|20|l . and the percolation criterion of Add. IdI 
Assuming that the conductivity saturates to Goo beyond 
the scale Vi ns , we may estimate Goo & G(V g = Vi ns ) from 
which we obtain the amplitude of the memory dip as 
AG = Goo - G(0) w G(V- ns ) - G(0). 

The smearing of the density of states due to new carri- 
ers has a minor effect on the non-equilibrium conductiv- 
ity since it mostly affects energy scales on the order of the 
temperature whereas the energy range probed by variable 
range hopping is much larger. However, the shift of the 
chemical potential, A/i is the crucial out-of-equilibrium 
feature which leads to the increase of the conductivity. 

The results we obtained from the percolation treat- 
ment confirm the general assertion 8 , that the conductiv- 
ity always increases with \V g \. In Fig. [S] we plot the rel- 
ative change of the conductivity AG(V^)/Goo for high 
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FIG. 7: (Color online) The relative amplitude of the conduc- 
tivity dip [G(Vins) — G(0)]/G(Vins) as a function of tempera- 
ture for various values of z. Inset: Low temperature behavior 
of the relative amplitude for the same values of z. 
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0.01 




FIG. 8: (Color online) Amplitude of the conductivity dip 
G(Vms) — G(0), as a function of temperature for z — 0.2, 
0.4 and 0.5. Less resistive films (larger z) exhibit a clear 
maximum at relatively high temperatures. 



temperatures T g > T > Ec- This is the temperature 
regime in which most experiments on indium-oxide films 
are performed. The cusp width decreases linearly with 
temperature while its amplitude increases, an asymptotic 
analysis yielding 



AG Ec_ 
~~G T ' 



(24) 



see Fig. □ 

At low temperature, T < Ec, the cusp width increases 
quadratically with T. In contrast to the high tempera- 
ture regime, the adiabatic percolation treatment com- 
bined with the instability criterion predict AG/Goo to 
increase with temperature as 



AG 
~~G~ 



T 



T, 



ES 



(25) 



The full functions AG(V g )/G OQ are shown in Fig. 

The absolute amplitude of the dip, AG, is found to be 
a non-monotonic function of temperature, as shown in 
Fig-El The non-monotonicity is more pronounced in less 
resistive films (see the curves for the localization param- 
eters z = 0.4, 0.5) where a clear peak appears in AG at 
Tmax> while for more resistive films (z = 0.2) this feature 
is hard to discern. Very similar non-monotonic behavior 
of AG was observed experimental!}^. 



V. DISCUSSION 

A. Low temperature behavior 

The non-monotonic behavior of the relative amplitude 
predicted by Eqs. (|24I25|I was not observed in experi- 
ments so far. Even though this might be due to the fact 



that most experiments are performed at high tempera- 
tures, we believe the prediction AG/G — > for T — ► 
to be an artifact for the following reasons: (i) One can 
trace back the origin of the prediction AG/G ~ T 1 ' 2 to 
the conservation of the total deficit of density of states 
under the application of a gate voltage. In our approxi- 
mation, this is a simple consequence of only shifting and 
convoluting the density of states. However, it is likely 
that for V ~ Vi ns the neglected fast relaxation processes 
destroy this exact conservation, which would lead to the 
saturation of AG/G at low temperatures, (it) The as- 
sumption that the asymptotic conductivity Goo is well 
estimated by the non-equilibrium conductivity at Vi ns is 
probably incorrect at very low temperatures. In order 
to illustrate this point, let us consider a large gate volt- 
age Vl with Vins < Vj, < V^„cm, such that the shift of 
the chemical potential is of the order of Ec- This will 
take the system into a high energy state where fast relax- 
ation processes lead to the formation of a new Coulomb 
gap. The sites within the new gap will be mostly dif- 
ferent from the sites in the old one, and it is likely that 
the total deficit of density of states in this new Coulomb 
gap is initially smaller than the one in the equilibrium 
state. More precisely, one may expect that the linear 
slope acgg of the density of states in the new configura- 
tion is slightly larger than that in the ground state, and 
the corresponding Efros-Shklovskii temperature Tes is 
smaller accordingly (see App. [DJ. This would suggest 
that AG/G scales like T -1 / 2 and thus increases with de- 
creasing temperature. In this temperature regime, one 
may expect an additional increase of conductivity with 
increasing gate voltage, even beyond the instability point 
Vi ns . The description of the conductivity in that regime 
would require to take into account partial relaxation pro- 
cesses which goes beyond the present approach. However, 
it remains open whether the VT behavior Eq. (|25|l ap- 



9 



plies in an intermediate range of temperatures. 



B. Comparison to experiments 

1. The memory dip 

As described in Section llV CI the adiabatic percolation 
approach reproduces well the temperature dependence of 
several key features of the memory dip. In particular, we 
showed that one may infer the width of the quasiparticle 
Coulomb gap from a careful study of the temperature 
dependence of the width. Moreover, from Eq. (|21|l we see 
that the width of the memory dip should be proportional 
to the bare density of states vq, which increases with 
carrier concentration as discussed in Section IIII Al In 
the impurity band regime, n c < Tlx, one expects v$ ~ 
r?J^ , while for large carrier concentration the density of 
states crosses over to a free electron-like behavior Uq ~ 

1/3 

n c ■ This scenario agrees rather well with the observed 
bending of the dip width as a function of carrier density 
in indium-oxide filmat 

At high temperatures, the percolation approach pre- 
dicts a decrease of the cusp amplitude like I/T, which is 
weaker than what is usually observed in experiments^*^. 
This difference may again originate from our neglect of 
spatial correlations in 5<j>i discussed above. Even more 
likely is the scenario that our assumption of homoge- 
neously glassy samples breaks down at higher temper- 
atures. Indeed, if only a few rare regions with stronger 
disorder remain glassy their effect on the out of equi- 
librium conductivity might be strongly reduced due to 
shortcuts by non-glassy regions. 

Let us briefly discuss the effect of varying the disorder 
strength or applying a magnetic field. Experimentally, 
both changes seem not to affect the width of the memory 
dip. Furthermore, even strong magnetic fields change 
the amplitude AG/G only slightly. These experimental 
observations are quite surprising since both the variation 
of disorder and magnetic field affect the conductivity G 
itself appreciably. 

The first observation can be naturally explained within 
our picture if we assume that changing the disorder (by 
annealing) or applying a magnetic field mostly affects the 
localization length without altering the bare density of 
states vq. Since the instability criterion originates from a 
static consideration within the classical model Hamilto- 
nian it is insensitive to the localization length. The 
width of the memory dip remains thus constant under 
variations of disorder or magnetic field. 

The near constancy of AG/G with magnetic field is 
more subtle. In terms of the percolation approach the 
variation of the localization length simply changes the 
parameter z, while leaving Eq fixed. From Fig. [3 it can 
be seen that the effect of z on AG/G is indeed relatively 
small, while the corresponding change of AG (cf. Fig. [SJ 
is much more important. 



2. Memory of temperature 

An interesting effect of temperature memory was re- 
ported both for the indium-oxide filmsi^, and for films 
of granular aluminum^. After equilibration at V g = 
and temperature To, the system is quenched to T\ < To, 
and the conductivity is probed as a function of gate volt- 
age before the sample equilibrates. In this protocol, the 
anomalous field effect maintains the characteristic width 
of the initial temperature To for a rather long time before 
narrowing down. 

This effect can be understood in terms of the instability 
criterion proposed above: The energy minima or valleys 
in which the electron glass typically settles at tempera- 
ture To will be stable upon injection of additional carri- 
ers up to the critical density nr - After a temperature 
quench the system will in general remain in this valley, 
and the stability threshold of the higher temperature will 
be preserved temporarily. Finally, slow relaxation pro- 
cesses allow the system to dig itself deeper down into a 
subvalley of the energy landscape, since the thermal fluc- 
tuations are now smaller. These lower lying states will 
have a reduced stability threshold, as will be reflected by 
a smaller dip width T. 

C. Observability of glassy effects in doped 
semiconductors 

So far, slow electronic relaxation was observed only 
in very few cases of moderately doped semiconductors 5 . 
The natural question thus arises as to why glassy effects, 
such as in indium-oxide, are not more frequently encoun- 
tered. The reason is most likely that many standard 
semiconductors have relatively low carrier concentrations 
with well localized electrons. Such systems are described 
by the classical impurity band model whose glass transi- 
tion in 3D is suppressed by a small numerical factor, see 
Eq. Estimating the Efros-Shklovskii conductivity 

at that temperature scale, one finds 

/ , P 2 \V2 / r (3D) -1/3 \ 1 ' 2 

toWRo) ,( cr _f_) , 

(26) 

which is very large even when the localization length ap- 
proaches the inter-impurity distance. This makes the 
detection of glassy effects in the hopping conductivity 
nearly impossible, due to the intrinsically large noise in 
such systems. However, glassiness should still be observ- 
able in static quantities, such as the capacitance mea- 
surements of Ref£. 

In amorphous semiconductors with relatively high car- 
rier concentration glassy effects as described in this pa- 
per should generally be observable. The same is true 
for doped semiconductors sufficiently close to the metal 
insulator transition. 

Once the localized wavefunctions start to overlap sig- 
nificantly (z > 1) one may expect the nature of the glass 
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phase to change and finally disappear completely. Many 
recent experiments probing the Coulomb gap are actu- 
ally carried out in this regime 4647 ' 59 . They reveal very 
interesting quantum critical behavior associated with the 
metal insulator transition, but have not thoroughly in- 
vestigated the glassy aspects of the samples so far. First 
attempts towards a theoretical description of glassiness 
in the regime close to the transition were undertaken in 
Refsiffii 6 ^. At this point it remains an interesting open 
question, both theoretically and experimentally, whether 
the onset of metastability and glassiness coincides with 
the transition to the insulator, and, if so, what role the 
glassy freezing plays in the physics of the metal-insulator 
transition. 



VI. CONCLUSION 

We have analyzed the memory effect in electron 
glasses. The non-equilibrium conductivity was calcu- 
lated within a percolation approach, assuming the lo- 
cal metastability of the glass state. This allowed us to 
describe the anomalous field effect quantitatively, repro- 
ducing many of the experimental characteristics observed 
in indium-oxides and granular aluminum. We have pro- 
vided a simple physical picture for the voltage scales at 
which the memory dip saturates and erasure of memory 
occurs, respectively. We argue that the saturation scale 
increases with temperature, its dependence on temper- 
ature reflecting the characteristics of the Coulomb gap. 
We have predicted the ratio of the two voltage scales as 
a function of temperature and carrier density, which can 
be tested in experiments. 
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APPENDIX A: THE BARE DENSITY OF STATES 
OF SEMICONDUCTORS 

In this appendix we derive approximate expressions for 
the bare density of states in semiconductors and granular 
metals. 

In the literature, two standard models for semiconduc- 
tors with localized electrons have been considered, see 
Refsi2ii22i for a review. The "classical impurity band" 
model refers to lightly doped, partially compensated 
semiconductors where all carriers are localized within 



a Bohr radius around majority impurities. Due to the 
Coulomb interactions with randomly distributed charged 
impurities, the on-site energies q of these localized states 
are scattered over a range of the order of the nearest 
neighbor interactions, e 2 / nr c , where r c = n c ^ 3 is the 
average distance between carriers, n c is the carrier den- 
sity (uncompensated dopant concentration) and k is the 
host dielectric constant. Accordingly, the bare density 
of localized states (which neglects Coulomb interactions 
between the localized carriers) is of the order of 

Vq = c — = — n 2 / 3 , (impurity band). (Al) 

e z jK,r c e z 

The second frequently considered model describes 
amorphous semiconductors in which the disorder of the 
on-site energies is due to strong local inhomogeneities. 
In this case, their scatter is usually much larger than that 
introduced by Coulomb interactions with impurities. As 
a consequence, the bare density of states is lower than 
(|A1|) . as schematically illustrated in Fig.|5| 

The set of localized states does not need to fill the 
whole region between the valence and the conduction 
band. In the case of indium-oxide (both amorphous 62 
and crystalline 6 ^), it has been established that the local- 
ized states form a tail joining the conduction band at the 
mobility edge. Furthermore, it was found that at suffi- 
ciently high carrier densities, the density of states in the 
range of localized states is in surprisingly good agreement 
with free-electron estimates, 

vq ~ 2/3 ' ( high densit y' n c > n x )- (A2) 

% n c /2m 

This reflects the fact that the kinetic energy Eki n — 
Ti 2 kp/2m ~ h 2 n'J' i jlm of the localized wavefunctions 
dominates over the effects of inhomogeneities in the elec- 
trostatic potential. Note, however, that a crossover to the 
regime of dominant Coulomb interactions (Eq. IjAlfl l is 

to be expected around n c ~ nx where E^in ~ e 2 n^ 3 /k, 
i.e., 

In granular metals, the role of the localized sites in 
the model is taken by the grains. As a consequence 
of impurities and the disorder in the size and the ar- 
rangement of the grains, the cost to introduce one more 
particle on a grain is a random quantity of the order of 
the typical charging energy Eq or the level spacing 5 
in the grain, whichever is larger. Usually the charging 
energy will dominate, unless the grains are very small 
or the effective dielectric constant of the metallic film is 
very large. The random on-site energies q entering the 
Hamiltonian Q are therefore scattered with a typical 
width W = max [Ec, S], and the effective bare density of 
states can be estimated as 

v ~ ft; — ft i (granular metals) (A4) 

max|fcc, o\ 
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FIG. 9: Schematic view of the bare density of states (neglect- 
ing electron-electron interactions) in different classes of semi- 
conductors: a) classical impurity band (lightly doped semi- 
conductors) b) strongly disordered amorphous semiconductor 
(e.g., amorphous germanium), c) doped, disordered semicon- 
ductor with localized band tails at the bottom of the conduc- 
tion band (e.g., indium-oxide). 



which is typically a few times smaller than the (2D) den- 
sity of states in a bulk metal. 



APPENDIX B: REDUCTION OF HIGH DENSITY 
SYSTEMS TO THE STANDARD MODEL 

Here we examine under which condition a high den- 
sity system can be described by the classical Hamilto- 
nian . The aim is to consider only a strip of localized 
states of width AE around the chemical potential, and 
to work with an effective model of occupied and empty 
levels within this strip. In this approximation, the carri- 
ers localized in states of lower energy are considered inert 
in the sense that they do not hop to other sites. Notice 
however, that such "core" electrons may still have fairly 
extended wavefunctions, and therefore contribute to the 
polarizability of the medium, renormalizing the host di- 
electric constant. 

The mapping of such an energy strip to a model of 
point-like localized states is consistent provided that 
(i) the states within the strip do not overlap spatially, 

AEv^ D < 1, (Bl) 

and that (ii) the typical variations Sip of the electrostatic 
energy due to rearrangements of particles within the strip 
do not exceed the width AE of the strip, 

2 

5<t> ~ {v AE) 1/D — < AE. (B2) 



T' 


i 


\ ln(R, 


~T-' M rp 


L 2D-3D X 


' L g 




ln(R)~T-" 3 \ 






/ /' ln(R)~T-" 2 






a 


cr f 



FIG. 10: (Color online) Phase diagram of electron glasses 
as a function of temperature and film thickness. The dash- 
dotted lines labeled Tx red (light grey) and T^d-ad purple 
(dark grey) indicate the crossover between different hopping 
regimes. The solid lines separate the ergodic high tempera- 
ture phase from the glass phase where memory effects and ag- 
ing are observable. The dashed line indicates the temperature 
and energy scale below which the density of states assumes a 
linear shape characteristic for two dimensions. 

The conditions I|B1|) and (|B2|) can be satisfied simultane- 
ously if 

z = -v Q t D - x < 1, (B3) 

K 

or, in other words, if the level spacing within a local- 
ization volume is larger than the Coulomb interaction 
strength on the scale of the localization length. 

APPENDIX C: THE CROSSOVER FROM 2D TO 

3D 

In this appendix we discuss the crossover thickness 
d cr below which a sample should be considered two- 
dimensional. In particular, we show that the crossover 
from a bulk sample to a film occurs around a thickness 

dcr ~ Ke 2 /«r 1/2 , (CI) 

both with respect to transport characteristics and the 
glass transition. [Using typical values for indium-oxide 
films (l>q w 10 32 er(? -1 cm -3 and k m 30^) one finds d cr ~ 
100A which is of the order of the typical film thickness 
(d = 50 — 200A) in most glassy experiments.] 

In a thin film, Mott's variable range hopping law Q 
crosses over from the 3-D form (log(i?) ~ T -1 / 4 ) to the 
2-D form (log(i?) ~ T -1 / 3 ) when the hopping length 
becomes of the order of the film thickness, which yields 
the crossover temperature 

T 2D - 3D ~ T§ D \t/d) A ~ Tf D) /{u d 2 e 2 ln)\ (C2) 

A subsequent crossover to the Efros-Shklovskii law takes 
place at T {2D) = T { ^ D) \v d 2 e 2 / '«) , where T^ D) denotes 
the Mott to Efros-Shklovskii crossover temperature for a 
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bulk sample. The intermediate regime with a 2D-Mott's 
law is observable only if d < d cr . 

The crossover from a 3D to a 2D glass transition occurs 
when the typical distance between thermally active sites 



at T = 
when 



Tg 3D ^ becomes equal to the film thickness, i.e., 



Rr a 



(C3) 



One can easily check that these expressions all become 
of the same order when vod 2 e 2 /k = (d/d cr ) 2 w 1. Notice 
that for films with d < d cr , the glass transition temper- 
ature @ decreases roughly linearly with thickness since 

Tg ~ fg 2 ~ dv{). The phase diagram as a function of 
temperature and film thickness is summarized in Fig. 1101 



APPENDIX D: PERCOLATION THEORY OF 
HOPPING CONDUCTIVITY 

In this appendix we review the percolation theory of 
hopping conductivity. We consider the network of Miller- 
Abrahams resistors formed by pairs of sites i and j. In 
the vicinity of a given low temperature metastable state 
of the electron glass, the effective resistance of this link 
is approximately given by 



i? exp(-2ry/^ + e ij /T) 



where 



f lEi-Ejl-e 2 /^, if Ei -Ej 
\ max {I^U^I}, KEi-Ej 



< 
> 



(Dl) 



(D2) 



and Ei is the energy (with respect to the chemical po- 
tential) to remove or add a particle at the site i in the 
particular metastable state at hand. More precisely, the 
energies Ei refer to the excitation of quasiparticles (or 
polarons^i) that carry the hopping current. 

In order to find the least resistive percolating path in 
the resistor network we follow the procedure proposed 
by Efros et ali^I: We consider only resistors with Rij < 
i?oexp(xc) to be active and associate to each of them a 
disk or ball with diameter . We finally determine the 
threshold value of Xc for which the set of disks percolates. 
The value Rq exp(% c ) is a good estimate of the resistivity 
to exponential accuracy. 

To solve this problem analytically, one needs to know 
the probability F(u>,r) per unit energy and volume to 



find a pair of sites (i,j) with = r and ej 



Un- 



der the assumption that the site energies Ei are inde- 
pendently distributed according to a single-quasiparticle 
density of states p(E), we obtain 



F(u,r) 



1 



p(E 1 )p(E 2 )S(e 12 -u)dE 1 dE 2 . (D3) 



With the help of the pair distribution function the above 
percolation problem reduces to that of a set of balls with 
different radii. Assuming that the critical volume fraction 
of the balls, Qd, is an approximate invariant of temper- 
ature, only dependent on dimensionality, we finally have 



to solve the equation 



J dojd D rV D (£} D F(u,r)9( X c - £ - j),(D4) 



where Vd is the volume of a £>-dimensional unit sphere. 
The above invariance principle yields 02 ~ 1.26 64 and 
O3 0.23>2»). Comparison with other percolation cri- 
teria, in particular in the Mott regime22iSiS& I indicate 
that in 2D a slightly smaller value 82 ~ 1 yields results 
closer to the numerically found percolation threshold. In 
the main part of the paper we therefore used the latter 
value. 

In order to efficiently implement the percolation crite- 
rion for an arbitrary density of states, as obtained, e.g., 
after the sudden application of a gate voltage, it is con- 
venient to introduce the functions 



F ph (E) = 2[ dep(e)p(e-E) 
Jo 

F pp (E)=p(E) f dep(e) 
Jo 



F hh {E) = p(-E) I dep{e) 

-E 



(D5) 
(D6) 
(D7) 



and Q a p(E) = L deF a p{e) with a,/3 £ {p, h}, in terms 
of which the percolation criterion (|D4|I can be rewritten 
as 



e D = 



'rV D ( - 



(D8) 



[% P {T( X c - 2r/0) + $hh(T( X c - 2r/0) 
+% h (T( Xc - 2r/0 + e 2 /nr) - <S> ph (e 2 / nr)] . 

For completeness we report the standard expressions 
that one obtains in the limiting case of a constant density 
of states (p(e) = vq) and high temperatures (T ^> Tx, 
Mott regime), and in presence of an Efros-Shklovskii 
pseudogap, p(e) = a D (n/ 'e 2 ) D t D ~ x , (T < T x )- In these 
cases the above criterion is readily evaluated analytically 
and yields the threshold values 



D s 1/(D+1) 



, (Mott) 



(D9) 



Xc = 



Ik 

T 

rpD \ !/ 2 

, (Efros-Shklovskii) (D10) 



with 



AD) _ 

M — 



r {D) 

voH D 



T 



•(D) 



ES 



c 



<(D) 
ES 



(Dll) 
(D12) 



The numerical constants effi depend on the value of 
D , while C^J increases with the value of Q£>/a D . In 
the main text we used <d 2 = 1- 



13 



1 J. Davies, P. Lee, and T. Rice, Phys. Rev. Lett. 49, 758 
(1982). 

2 M. Griinewald, B. Pohlmann, L. Schweitzer, and D. Wiirtz, 
J. Phys. C 15, L1153 (1982). 

3 J. Davies, P. Lee, and T. Rice, Phys. Rev. B 29, 4260 
(1984). 

4 M. Pollak, Philos. Mag. B 50, 265 (1984). 

5 D. Monroe, A.C. Gossard, J.H. English, B. Golding, and 
W.H. Haemmerle, Phys. Rev. Lett. 59, 1148 (1987). 

6 M. Ben-Chorin, Z. Ovadyahu, and M. Pollak, Phys. Rev. 
B 48, 15025 (1993). 

7 A. Vaknin, Z. Ovadyahu, and M. Pollak, Phys. Rev. Lett. 
81, 669 (1998). 

8 A. Vaknin, Z. Ovadyahu, and M. Pollak, Phys. Rev. Lett. 
84, 3402 (2000). 

9 V. Orlyanchik and Z. Ovadyahu, Phys. Rev. Lett. 92, 
066801 (2004). 

10 Z. Ovadyahu, Philos. Mag. B-Phys. Condens. Matter Stat. 
Mech. Electron. Opt. Magn. Prop. 81, 1225 (2001). 

11 Z. Ovadyahu and M. Pollak, Phys. Rev. B 68, 184204 

(2003) . 

12 A. Vaknin, Z. Ovadyahu, and M. Pollak, Phys. Rev. B 65, 
134208 (2002). 

13 C.J. Adkins, J.D. Benjamin, J.M.D. Thomas, J.W. Gard- 
ner, and A.J. McGeown, J. Phys. C 17, 4633 (1984). 

14 E. Bielejec and W. Wu, Phys. Rev. Lett. 87, 256601 (2001). 

15 T. Grenet, Eur. Phys. J. B 32, 275 (2003). 

16 G. Martinez-Arizala, D.E. Grupp, C. Christiansen, A.M. 
Mack, N. Markovic, Y. Seguchi, and A.M. Goldman, Phys. 
Rev. Lett. 78, 1130 (1997). 

17 G. Martinez-Arizala, C. Christiansen, D.E. Grupp, N. 
Markovic, A.M. Mack, and A.M. Goldman, Phys. Rev. B 
57, R670 (1998). 

18 T. Grenet, phys. stat. sol. (c) 1, 9 (2004). 

19 A. Vaknin, Z. Ovadyahu, and M. Pollak, Europhys. Lett. 
42, 307 (1998). 

20 S. Bogdanovich and D. Popovic, Phys. Rev. Lett. 88, 
236401 (2002). 

21 J. Jaroszynski, D. Popovic, and T. M. Klapwijk, Phys. 
Rev. Lett. 89, 276401 (2002). 

22 M. Pollak, Disc. Faraday Soc. 50, 13 (1970). 

23 A. L. Efros and B. I. Shklovskii, J. Phys. C 8, L49 (1974). 

24 A. L. Efros, J. Phys. C 9, 2021 (1976). 

25 A. Pastor and V. Dobrosavljevic, Phys. Rev. Lett. 83, 4642 
(1999). 

26 M. Miiller and L. B. Ioffe, Phys. Rev. Lett. 93, 256403 

(2004) . 

27 S. Pankov and V. Dobrosavljevic, Phys. Rev. Lett. 94, 
046402 (2005). 

28 A. Zabrodskii, Philos. Mag. B-Phys. Condens. Matter Stat. 
Mech. Electron. Opt. Magn. Prop. 81, 1131 (2001). 

29 Physical aging in Amorphous Polymers and other materi- 
als, edited by L. C. E. Struick (Elsevier, Houston, 1978). 

30 P. NordBlad and P. Svedlindh, in Spin Glasses and Ran- 
dom Fields, edited by A. P. Young (World Scientific, Sin- 
gapore, 1998). 

31 A. Efros and B. Shklovskii, in ElectronElectron interaction 
in disordered systems, edited by A. Efros and M. Pollak 
(North-Holland, Amsterdam, 1985). 

32 B. Shklovskii and A. L. Efros, Electronic properties of < 
semiconductors (Springer, Heidelberg, 1984). 



"Polaron" excitations consist in the change of occupation 
on one site, and a small number of short hops of neigh- 
boring particles that decrease the energy cost for the local 
charge excitation—. 

34 V. Ambegaokar, B. I. Halperin, and J. S. Langer, Phys. 
Rev. B 4, 2612 (1971). 

35 B. I. Shklovskii and A. L. Efros, Zh. Eksp. Theor. Fiz. 60, 
867 (1971). (Sov. Phys. JETP 33, 468 (1971)). 

36 M. Pollak, J. Non-Cryst. Solids 11, 1 (1972). 

37 A. L. Efros, V. L. Nguyen, and B. I. Shklovskii, Solid State 
Commun. 32, 851 (1979). 

38 A. Aharony, Y. Zhang, and M. P. Sarachik, Phys. Rev. 
Lett. 68, 3900 (1992). 

39 Y. Meir, Phys. Rev. Lett. 77, 5265 (1996). 

40 J. Massey and M. Lee, Phys. Rev. Lett. 75, 4266 (1995). 

41 M. Lee, J. Massey, V. Nguyen, and B. Shklovskii, Phys. 
Rev. B 60, 1582 (1999). 

42 W. Teizer, F. Hellman, and R. Dynes, Phys. Rev. Lett. 85, 
848 (2000). 

43 B. Sandow, K. Gloos, R. Rentzsch, A.N. Ionov, and W. 
Schirmacher , Phys. Rev. Lett. 86, 1845 (2001). 

44 V. Butko, J. DiTusa, and P. Adams, Phys. Rev. Lett. 84, 
1543 (2000). 

45 V. Y. Butko, J. F. DiTusa, and P. W. Adams, Phys. Rev. 
Lett. 85, 162 (2000). 

46 V. Y. Butko and P. Adams, Nature 409, 161 (2001). 

47 M. Lee, Phys. Rev. Lett. 93, 256401 (2004). 

48 F. Ladieu, F. Bert, V. Dupuis, E. Vincent, and J. Ham- 
mann, J. Phys. C 16, S735 (2004). 

49 A. Perez- Garrido, M. Ortuno, A. Diaz-Sanchez, and E. 
Cuevas, Phys. Rev. B 59, 5328 (1999). 

50 D. Menashe, O. Biham, B. Laikhtman, and A. Efros, Eu- 
rophys. Lett. 52, 94 (2000). 

51 D. Menashe, O. Biham, B. Laikhtman, and A. Efros, Phys. 
Rev. B 64, 115209 (2001). 

52 D. Tsigankov, E. Pazy, B. Laikhtman, and A. Efros, Phys. 
Rev. B 68, 184205 (2003). 

53 D. Grempel, Europhys. Lett. 66, 854 (2004). 

54 A. B. Kolton, D. R. Grempel, and D. Dominguez, Phys. 
Rev. B 71, 024206 (2005). 

55 E. R. Grannan and C. C. Yu, Phys. Rev. Lett. 71, 3335 

(1993) . 

56 M. Overlin, L. Wong, and C. Yu, Phys. Rev. B 70, 214203 
(2004). 

57 M. Miiller and E. Lebanon, cond-mat/0508703 to be pub- 
lished in J. Phys. France IV (2005). 

58 Z. Ovadyahu, private communication. 

59 E. Helgren, N. Armitage, and G. Gruner, Phys. Rev. B 69, 
014201 (2004). 

60 D. Dalidovich and V. Dobrosavljevic, Phys. Rev. B 66, 
081107 (2002). 

61 V. Dobrosavljevic, D. Tanaskovic, and A. Pastor, Phys. 
Rev. Lett. 90, 016402 (2003). 

62 Z. Ovadyahu, Phys. Rev. B 47, 6161 (1993). 

63 O. Cohen and Z. Ovadyahu, Phys. Rev. B 50, 10442 

(1994) . 

64 V. L. Nguyen, Sov. Phys. Semicond. 18, 207 (1984). 

65 C. H. Seager and G. E. Pike, Phys. Rev. B 10, 1435 (1974). 

66 G. E. Pike and C. H. Seager, Phys. Rev. B 10, 1421 (1974). 



